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ABSTRACT 

It is strongly believed that Andromeda's double nucleus signals a disk of stars re- 
volving around its central super-massive black hole on eccentric Keplerian orbits with 
nearly aligned apsides. A self-consistent stellar dynamical origin for such apparently 
long-lived alignment has so far been lacking, with indications that cluster self-gravity 
is capable of sustaining such lopsided configurations if and when stimulated by ex- 
ternal perturbations. Here, we present results of N-body simulations which show un- 
stable counter-rotating stellar clusters around super-massive black holes saturating 
into uniformly precessing lopsided nuclei. The double nucleus in our featured experi- 
ment decomposes naturally into a thick eccentric disk of apo-apse aligned stars which 
is embedded in a lighter triaxial cluster. The eccentric disk reproduces key features 
of Keplerian disk models of Andromeda's double nucleus; the triaxial cluster has a 
distinctive kinematic signature which is evident in HST observations of Andromeda's 
double nucleus, and has been difficult to reproduce with Keplerian disks alone. Our 
simulations demonstrate how the combination of eccentric disk and triaxial cluster 
arises naturally when a star cluster accreted over a pre-existing and counter-rotating 
disk of stars, drives disk and cluster into a mutually destabilizing dance. Such accre- 
tion events are inherent to standard galaxy formation scenarios. They are here shown 
to double stellar black hole nuclei as they feed them. 
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1 INTRODUCTION 

Over the years, the nucleus of the Andromeda galaxy (M31) 
went from being asymmetric, to doubling, then tripling. The 
asymmetry was first noted in the balloon-born Stratoscope 
observatory l|Light et al.lll974l ). It was photometrically re- 
solved into a double nucleus by Hubb le Space Telescope 
(HST) observations |Lauer et alj|l993j ): asymmetry hence 
turned into lopsidedness, with a luminous feature (referred 
to as PI) shining a few parsecs away from a dimmer one 
(referred to as P2), which is closer to the center of the host 
bulge. The resolved double nucleus had already been sus- 
pected to host a Super-Massive Black Hole (SMBH), located 
somewhere close to P2, with a neighborhood that was known 
to sh ine brighter than the rest of the nucleus in t he ultravi- 
olet (|Dressler fc Richstondll988l : H ormcn dvl ll988l ). Detailed 
HST spectroscopy added a third component (known as P3), 
a disk of young massive stars, whose size and rotation speed 
rules out viable alte rnatives to a centra l SMBH of ~ 10 8 
solar masses (M® ) (|Bender et al.ll2005l ). 

In the c urrently favored model of M31's lopsided dou- 
ble nucleus (|Tremaind Il995l : IPeiris fc Tremainel l2003h . the 
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brighter peak PI (|Lauer et al.l Il99&t ) is thought to coin- 
cide with the common apo-centric region of an eccentric 
disc of stars revolving on nearly apse-aligned Keplerian el- 
lipses, around M31's central super-massive black hole (here- 
after SMBH). Similar eccentric discs are thought to under- 
lie lopsided nuclei detected around SMBH s in the centers 
of ne arby galaxies l|Lauer et al.l 1 199*61 , 120051 ; iGultekin et al.l 
l201ll) . That such kinematic configurations can be sustained 
in self-gravitating disks around SMBHs was clarified in a 
series of dynamical studies. Indeed, investigation of modes 
of hot nearly-Keplerian disks indicated that slow (m=l , 
lopsided) modes can be stably excited l|Tremainel 1200 ll ). 
This conclusion was corroborated by results of planar N- 
body simulations of disks around SMBHs which showed 
that, when started in asymmetric conditions, disks can 
relax into l ong-lived un i forml y precessing lopsided con - 
figurations ([Bacon et al.l l200ll ; Ijacobs fc Sellwoodl l200ll ). 
Subsequently, razor-thin, self-gravitating, lopsided equilib- 
ria were constructed, yielding encouraging agreement with 
photometric and kinematic observations of M31's nucleus 
l|Salow fc Statlerl l200ll ; ISambhus fc Sridharl |2002| ). These 
and related studies left open questions about how such 
modes are excited, how they ultimately saturate into global 
lopsided configurations, and whether constructed equilibria 
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were stable or not. The early suggestion (|Tremainei 1 19951 ') 
that an initially circular disk could become eccentric under 
the influence of dynamical friction with the host bulge has 
not been thoroughly explore d. A recently proposed scenario 
l|Hopkins fc QuataertlfcoiOal lbl) has eccentric disks forming 
in th e notoriously complex -and poorly understood fjSilld 
2011)- environment which couples star formation and its 
feedback to gas accretion in the early stages of SMBH 
growth. 

In this work, we opt for a minimalist and dynamically 
self-consistent route to three-dimensional, lopsided nuclei, 
out of counter-rotating (CR hereafter) collisionless stellar 
distributions dominated by SMBHs. Such CR distributions 
are prone to violent m= l instabilities in the p r esence of mod - 
erate counter-rotation (| Jog fc Combesll2009l ; lToumal f2002). 
hence the suggestion that the lopsided structure in An- 
dromeda's nucleus may have been triggered by the accre- 
tion of a retrograde globular clust er in a pre-existing disk 
of stars (jSambhus fc Srid har 2002). With the intention of 
exploring the outcomes of such accretion events, we per- 
formed a series of N-body simulations of CR stellar distri- 
butions evolving in the sphere of influence of an SMBH. 
Below, we report on results which show unstable CR distri- 
butions evolving into stable, uniformly precessing lopsided 
nuclei. We follow the instability in our featured simulation 
from multiple complementary perspectives. We then probe 
the ensuing lopsided nucleus, and show that in addition to 
displaying all the key observational signatures which Ke- 
plerian disks are meant to model, it recovers asymmetries 
in the tail of Adrom eda's line of sight velocity distribution 
l|Bender et al.ll2005|)_ which a thick lopsided disk alone is un- 
able to reproduce! Peiris fc Trema inc 2003). 



2 COUNTER- ROTATING INSTABILITY: 
FROM LINEAR GROWTH TO 
SATURATION 

Our N-body simula tions were per formed with the parallel 
version of Gadget- 2 (Springcl 2005), an oct-tree code, which 
is popular with the cosmological structure formation com- 
munity. We pushed this tool to the limit of extreme mass 
ratios, one for which it was not specifically designed, one 
in which it did the desired job, albeit with stringent force 
accuracy, and time stepping criteria. With M31 in mind, all 
experiments have an SMBH with a mass of M. = 10 8 M© and 
a main disk compone nt (prograde in our convention) with a 
tenth of SMBH mass (jBender et al.ll2005l ; |Peiris fc Tremaind 
2003). An exhaustive exploration, with sufficient realism, of 
a scenario m which a CR cluster (10 5 - 1 6 M Q ) decays 
unde r dynamical friction to nuclear regions (iTremaine et alj 



1975), then gets disrupted (|Kim fc Morrial2003l ; iFuiii et alj 
20101 ) as it couples to the massive nuclear stellar disk, would 
have been forbidding to pursue in detail with our current 
computational resources. Instead we opted for an extensive 
exploratio n of a realist i c configuration which is known to be 
unstable (|Toum a 2002; S ridhar fc Sainill2010l ). and in which 
the disrupted CR cluster is modeled by overlaying the pro- 
grade disk with a retrograde disk of 1/10 its mass. Details 
on initial conditions for this experiment, along with Gadget- 
2 simulation parameters, are described in Appendix[X] vari- 
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Figure 1. The eccentricity-inclination dynamics of the particles 
up to 1 Myr is shown. In row a, we follow a top view of the 
prograde particles, from initial axisymmetric annular configura- 
tion, through a phase of growth of lopsidedncss, with growing 
mean eccentricity (0.41 around 1 Myrs), and alignment of ap- 
sides. Row b shows a top view of retrograde particles over the 
same period; they experience greater increase in mean eccentric- 
ity (around 0.73 by 0.72 Myrs), and show signs of apse alignment 
by 0.5 Myrs, before dispersing (by 1 Myrs) into a disky blob. A 
side view of retrograde particles in row c shows how the in-plane 
dispersal seen in row b reflects in projection a burst of out-plane 
dynamics which puffs the initially retrograde disk into a triaxial 
cluster of stars. The eccentricity-inclination dynamics of the full 
cluster is captured in row d, with prograde in black, and retro- 
grade in blue; the mean eccentricity of both populations grows as 
the instability develops; by 1 Myrs, the mean inclination of the 
prograde particles has hardly changed, while the highly eccentric 
retrograde population, is widely dispersed in inclination. 



ations on the fiducial experiment are discussed in Appendix 



B 



2.1 Dissecting the Instability 

CR configurations are known to be (linearly) unstable 
ijToum a 2002; S ridhar fc SaintolOh . We are here concerned 
with the (non-linear) saturation of this instability and shall 
first map its unfolding with complementary views of par- 
ticle dynamics (Fig. Q]). A top view of stars, in both pro- 
grade and retrograde populations (hereafter PP and RP 
respectively), shows them developing lopsidedness by cou- 
pling growing eccentricity to apse-alignment. The more mas- 
sive PP maintains its coherent, apse-aligned precession (Fig. 
03(a)) whereas the lighter RP dissolves gradually into a 
disky blob after 0.5 Myrs (Fig. [U(b)). Viewed from the 
side, the gradual in-plane dispersal of the RP occurs along 
with a dramatic excitation of out-of-plane motion; by 1 Myrs 
the RP unfolds into a triaxial structure (Fig. [U(c)). Scatter 
plots of particle eccentricities and inclinations (Fig. HJ(d)) 
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Figure 2. The linearly unstable, counter-rotating configuration 
saturates into a lopsided uniformly precessing disk. Panels a and 
b show the mean eccentricity <e> and inclination <I> of PP 
(black) and RP (blue). <e> grows, at near constant <I>, till 
about 0.5 Myrs, at which point <I> takes off, with both <e> 
and <I> peaking around 0.72 Myrs. Cycles of increasingly smaller 
amplitude follow this initial growth phase, leading to a saturated 
eccentric prograde disk which is shrouded by a triaxial halo. In 
panel c, the power of the four dominant modes in projected den- 
sity is displayed; an initially axisymmetric configuration (m=0) 
lends way to a dominant m=l mode, along with non-negligible 
m=2 and m=3 contributions. Triaxial dispersal of the RP is re- 
sponsible for the slight growth of the m=0 amplitude, and the re- 
lated slight decrease in the m=l mode amplitude after 0.72 Myrs. 
In panel d, we follow the precession rate of the m=l mode in time: 
initial large amplitude oscillations reflect back and forth libration 
of prograde and retrograde m=l excitations about the uniform 
precession state; these oscillations die out with the dispersal of 
the retrograde bunch, leaving an m=l mode which precesses at a 
near constant rate of 19 km s — 1 pc^ 1 . 



reveal how, by 0.72 Myrs, the now highly eccentric RP is 
well on its way to complete triaxial disruption. Tucked in 
this emerging cluster, the PP, which has absorbed much of 
the RP's (negative) angular momentum, heats up slightly 
in inclination as it consolidates its eccentric lump, then ad- 
justs slowly to whatever little angular momentum is left in 
the RP. 

A look at the mean eccentricity and inclination of both 
populations (Fig. [2] a, b) confirms these observations, as 
it reveals three distinct phases: Phase-I of near coplanar 
growth of the mean eccentricity of both populations which 
lasts till about 0.5 Myrs; Phase-II of continued mean eccen- 
tricity growth, but now coupled to growth in the mean in- 
clination of the RP, both of which reach a maximum around 
0.72 Myr; Phase-Ill, with cycles of increasingly small am- 
plitude leading to a near-steady regime. As emphasized in 
the caption of Fig. [2] transitions in mean eccentricity and in- 
clination are neatly imprinted on the modal structure of the 
projected density and the pattern speed of the dominant 
m=l mode, as the full cluster transitions from a thin ax- 
isymmetric initial state, to a saturated, lopsided, and slowly 
precessing mode (Fig. [21 c, d, and Fig. [3J). 

Saturation of the m=l mode is evidently correlated 
with the dispersal of the initially retrograde population into 
a triaxial cluster of high eccentricity orbits. The disper- 
sal appears to be associated with collective dynamics along 
eccentricity-inclination cycles. These cycles are excited when 
a population of stars (the retrograde population) develops 
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Figure 3. The projected density of the saturated configuration 
at the end of the simulation is followed over a precession cycle; 
a Pl-like peak is evident around apse-alignment; the distribution 
is practically invariant as it revolves uniformly over its 0.32 Myr 
cycle. 



sufficient eccentricity (through a CR instability) in the pres- 
ence of a dominant eccentric perturber (the more massive ec - 
centric and precessing prog rade disk) iTouma et ail (2009). 
One can presumably describe the incept ion of these cycle s 
in a generalized Kozai-Lidov framework (|Naoz et al.ll201ll ). 
but then one is left with the harder task of accounting for 
the dispersal of the population on these cycles. We gained 
valuable insight into dispersal then saturation by modeling 
a closely related process in a 2D analog of the 3D cluster 
(Kazandjian & Touma, in preparation). The model in ques- 
tion permits a blow by blow account of the mutual sculpting 
of planar prograde and retrograde populations in terms of 
capture in, t hen escape from a drifti ng, trapping region in 
phase-space l|Sridhar fc Toumal [l997h . It can naturally ex- 
plain the near-coplanar growth and apse alignment in phase 
I of 3D dynamics. Furthermore, it clarifies the collisionless 
damping process at work in the decaying cycles of phase III. 
How the eccentricity growth of phase I paves the way for in- 
clination growth in phase II cannot be captured in a planar 
analogue and shall await a more sophisticated treatment of 
the 3D deployment of unstable CR clusters. 



3 CONFRONTATION WITH M31: 
PRELIMINARY RESULTS 

The saturated distribution of stars, which is followed 
in projection, and over a full precessional cycle in Fig. 
[3] generalizes Keplerian disk models of M31's nucleus 
IPeiris fc TremainefeOOSt ) in two obvious ways. On one hand, 
the eccentric disk in our nucleus is shrouded by a lighter 
triaxial cluster, whereas currently preferred models tend to 
work with the eccentric disk alone: as discussed below, the 
additional triaxial cluster contributes a crucial improvement 
to model kine matics. On the other han d, successful eccentric 
disk models (|Peiris fc Tremaine] I2003I ) are essentially kine- 
matic in nature (their origin, and dynamical evolution re- 
maining uncertain), whereas ours is a dynamically stable 
nucleus to which an unstable configuration saturates. Given 
that the eccentric disk component in our saturated lopsided 
nucleus is ten times more massive than the triaxial cluster, 
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Figure 4. Photometry and kinematics of saturated configura- 
tions are probed, and a remarkable agreement with the kinemat- 
ics of Andromeda's nucleus is demonstrated. In the top row, 
we show the results of photometric and kinematic observations 
of our fiducial simulation: a double-peak is evident in the surface 
brightness distributiontop [left], with a brighter (Pl-like) peak on 
the left of the SMBH at the origin, and a fainter (P2-like) peak 
on its right; the zero point in the rotation curve top [center] 
is shifted towards brighter peak and slower rotation; the peak 
velocity dispersion is displaced towards fainter peak and faster 
rotation [right] ; rotation and LOSVD curves are shown without 
(dashed) and with (solid) contributions from the dispersed CR 
disk, which contributions introduce a signifi cant asymmetry in 
the ta ils. The best fit (non-aligned) models of lPeiris fc Tremainj 
(2003) (labeled PT03) with no bulge are also displayed as dot- 
dashed lines. A Quantitative agreement is demonstrated in the 
bottom row, center and right panels, between a dynamically 
similar copy of the featured simulation (obser ved under HST con - 
ditions) (solid line) and M31's kinematics l|Bender et al.|[2005f l 
(blue and red dots). 



the photometry and kinematics of this nucleus will naturally 
display all the significant qualitative features which eccentric 
disk models seek to explain (Fig. [4] top row): a double peak 
in surface brightness, an asymmetry in the rotation velocity 
(with an off-center zero velocity point), and a n off-centered 
peak in the line of sight velo city dispersion (| Bender et al.l 
l2005l : iPeiris fe Tremaine!l2003l ) (hereafter LOSVD). Pushing 
our luck, we seek a preliminary quantitative confrontation 
of simulation results with observations of M31's nucleus. In 
so doing, we work with a rescaled and dynamically sim- 
ilar version of our fiducial experiment (roughly speaking, 
we rescale from fiducial simulation length scale ~ 1 pc, to 
M31 nucleus scale, ~ 3 pc). The result ing Model Nucleus 
is then observed under H ST conditions (|Bender et al.llioOol : 
IPeiris fc Tremainell2003h (details in Appendix[£}. The out- 
come of this procedure is then overlaid with M31's kinemat- 
ics (Fig. |4l bottom row). We closely match the observed shift 
in the zero-point of rotation curve (x=-0.485 co mpared to 
x=-0.448 pc in observations (|Bender et al-l lioOS) ), the dif- 
ference in the peak rotation velocities (v=585 km s" 1 com- 
pared to v=532 ± 41 km s" 1 ), the magnitude of and shift 
in the peak in the LOSVD curve (x=+0.187, cr=+355.7 km 
s" 1 compared to x=0.298 pc, cr=+373 ± 48 km s" 1 ). 

The triaxial cluster of our Model endows its kine- 
matics with the marked asymmetry between the tails 
of the LOSVD, which is clearly observed in HST data 
|Kormendv fc Benderlll999l ; iBender et al.ll2005l ). and which 



has proven difficult to r ecover with eccentric disks alone 
jPeiris fc Tremaindliooll ') (compare PT03 and the "no CR" 
curves to the solid curve in the rightmost column of Fig. 
3J. Recall that this cluster consists of eccentric orbits which 
result from the triaxial dispersal of the initially retrograde 
cluster (rows b and d of Fig. [TJ). When superposed over the 
more massive lopsided disk, stars in this cluster will leave 
their kinematic signature mostly in the outer parts of the nu- 
cleus, which is where they linger the most on their eccentric 
orbits. That signature is particularly pronounced where the 
lopsided prograde disk contributes least, i.e. on the anti-Pi 
side of the nucleus, hence the resulting asymmetry. 

The initial ring-like configuration in our experiment is 
clearly deficient in stars in its central region. Hence, we did 
not expect for the resulting double nucleus to closely approx- 
imate the observed surface brightness profile in that region. 
This said, we were pleased to learn that a surface bright- 
ness profile, along a slit extending from brighter to fainter 
peak, has maxima around 13.3 and 14.2 mag.arcsec" 2 re- 
spectively, which compare favorably with the 13.2 and 13.57 
mag.arcsec -2 , observed by HST. 

Last but not least, we note that for the adopted similar- 
ity transformatiorQ the pattern speed of the lopsided mode 
scales down from the 19 km s -1 pc" 1 of the fiducial nucleus 
(Fig. H d), to the slower 4.3 km s" 1 pc" 1 of the Model 
Nucleus. At this relatively slow pattern speed, gravitational 
perturbations by the Model Nucleus are expe cted to drive 
and c onfine gas in a disk close to the SMBH (I Chang et alj 
l2007t ). When fueled by stellar mass loss from the nucleus, 
such a disk can regenerate starbursts at a rate which is con- 
sistent with the observed c ompact cluster of early type stars 
(P3) around M31's SMBH l|Bender et al.ll2005l ; Ichang et all 
2007). Thus, M31's P3 could very well be a direct conse- 
quence of PI and P2 in our dynamical Model. 



4 DISCUSSION 

Andromeda's kinematics shows a marked asymmetry be- 
tween the tails of its LOSVD which is clearly unaccounted 
for in Keplerian disk models, and appears to be cleanly re- 
solved by combining an eccentric disk with a triaxial cluster 
of highly eccentric orbits. To be sure, thorough modeling 
is called for before we can rigorously quantify the improve- 
ments of a disk-cl uster combination over t he thick eccen- 
tric disk model of lPeiris fc Tremainel <|2003r ). Still, we find 
it remarkable that this combination achieves, with minimal 
probing, the level of qualitative and quantitative agreement 
described above, and to see it emerge as a natural end prod- 
uct of our process ca n only strengthen t he case for CR stim- 
uli of double nuclei (|Lauer et al.l [2005). That case is made 
stronger when one learns that a rich variety of CR configura- 
tions are equally prone to developing lopsidedness (details in 
Appendix [B]). Experiments with CR point masses suggest 
that a cluster, initially on a circular trajectory, can drive 
a coplanar disk unstable provided its orbital radius is less 
than a critical radius (which is roughly equal to twice the 



In this transformation, time is rescaled by si , whenever length 
is rescaled by s. The adopted s ~ 2.7 accounts for the slower 
precession rate of the Model Nucleus. 
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mean radius of the disk); clusters on highly inclined circular 
trajectories can still drive a CR disk into a lopsided state; 
on the other hand, clusters with too large an initial eccen- 
tricity will end up with little negative angular momentum to 
drive a coplanar disk lopsided. The CR perturbation envis- 
aged in this work is delivered by a CR cluster which spirals 
deep into the nucleus b y dynamical friction with bulge stars 
(Tr emaine et al.|[l975h before it disrupts in the combine d 
tidal field of bulge and SMBH (|Quillen fc Hubbard! I2003T ) . 
The orbital radius at which the cluster disrupts depends crit- 
ically on the cluster's core density. For a cluster migrating 
on a near-circular trajector y into M31's nucleus, we estimate 
jQuillen fc Hubbard! [20o3 ) that a core density > 5 x 10 s 
Mq pc -3 is sufficient for the cluster to cross past the critical 
radius for instability. Mass s egregation and /or intermediate- 
mass black hole formation (|Kim Sz Morrisll2003l ; iFuiii et all 
can amplify core density significantly, thus further de- 
laying disruption till the cluster migrates into the overlap- 
ping configurations explored here. Such configurations may 
also arise when a cluster penetrates the nucleus on an ec- 
centric orbit, and disrupts upon close encounter with the 
SMBH. CR clusters will likely approach a pre-existing nu- 
clear disk on eccentric, and inclined trajectories; they may 
thus find themselves on the eccentricity-inclination cycles 
observed above, then disperse as they drive the disk in a 
saturated lopsided configuration. Given on one hand the 
strong l ikelihood of counter-r otating excitations in galactic 
centers (|jog fc CombesH2009r ). and on the other the robust- 
ness and efficiency of the proposed mechanism, double nuclei 
are likely to prove ubiquitous in stellar clusters dominated 
by super-massive black holes (|Lauer et alj|2005h . More gen- 
erally (and irrespective of whether or not a given lopsided 
nucleus results from a CR instability), the proposed mecha- 
nism can be deployed to customize triaxial equilibrium con- 
figurations with whic h to model obser ved kinematics of stel- 
lar black hole nuclei l|Alexanderll2005h (the Milky Way's in- 
cluded) and improve esti mates of the mass and feeding raten 
of th e black hole within (|Magorrian et al". Hl998l : lKormendvl 
I2004T ). Whether exploring realistic CR scenarios or tailor- 
ing triaxial equilibrium models, extensive numerical simu- 
lations of CR stellar systems with s tate of the art solvers 
|Fuiii et al.ll20ld ; iTouma et alj|2009h will surely contribute 
invaluable insights into the dynamics and structure of stellar 
black hole nuclei. 
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A METHODS 

We performed our si mulations usin g; the cosmological nu- 
merical tool Gadget-2 l|Springell2 005). Forces were evaluated 
using the oct-tree; cosmological, SPH, and grid capabilities 
were switched off. Although this narrowed down the param- 
eter space significantly, one had to be careful fine tuning 
parameters which control tree walk and time stepping. The 
depth of the tree to be explored is controlled by Gadget-2's 
parameter ErrTolForceAcc: daughter cells are opened until 
a certain criterion, set by ErrTolForceAcc, is satisfied, thus 
terminating the tree walk (we refer the reader to Gadget- 
2's user guide for more details about this criterion). Once 
the accelerations of all the particles are at our disposal, the 
positions of these particles are advanced with a time-step 
At = W2rjb/\a\, where |a|, r\ and e are the magnitude of the 
computed acceleration, ErrTolIntAccuracy (which controls 
the step-size) and the gravitational softening respectively. 
Setting ErrTolForceAcc to 0.001 and ErrTolIntAccuracy 
to 0.001 assured accurate enough simulations for our pur- 
poses. 

Our featured simulations consist of a massive prograde 
disk, on which a light retrograde disk is overlayed, with both 
residing in the sphere of influence of a central SMBH. Such a 
configuration is a likely outcome of the disrup tion of clusters 
that can migrate deep in the galactic c enters (|Gerhardll200ll : 
iKim fc Morris! l2003l : lFuiii et al.ll2010h . Initial positions and 
velocities of disk particles are picked from radial and vertical 
distri butions which are standard to disk dynamics simula- 
tions (jEmsellem fc Combesl H997h. Since Gadget-2 has diffi- 
culties in dealing with high mass density contrast, we further 
chose our disk to be annular in-order to avoid close encoun- 
ters with the SMBH. Thus, we introduced an inner and outer 
radial cut-off when sampling the disks. The SMBH is mod- 
eled as a Plummer-softened particle with softening length 
of b=0.01 pc. To minimize sources of instability which are 
not related to counter-rotation, disk components were sep- 
arately virialized before coupling them to each other or to 
perturbers. In the course of virializing, the more massive 
disks (prograde in our simulations) underwent short-lived 
gravitational instabilities, before relaxing to a hotter near- 
axisymmetric equilibrium state. 

The central SMBH with a mass of 10 s Mq is ten times 
more massive that the pre-existing disk of stars, itself ten 
times more massive than the disturbing CR disk, the mean 
eccentricities of the virialized prograde and retrograde disk 
were <e p > — 0.22 and <e r > = 0.25 respectively, with cor- 
responding standard deviations a p = 0.12 and cr r = 0.14. In 
both disks, the mean of the semi-major axis distribution was 
<a>~1.06 pc with a standard deviation of ~0.3 pc whereas 
the spread in inclination was around 8° . 

Simulation particles revolve around the central SMBH 
on perturbed Keplerian ellipses; associated osculating or- 
bital elements [namely semi- major axis (a), eccentricity (e), 
inclination (I), argument of periapse (g), and longitude of 
the node (h)] are simply recovered from particle position 
and velocity. In addition to particle orbit elements, our dy- 
namical analysis required the time evolving Fourier modes 



of the surface (projected) density of the CR cluster, and as- 
sociated potentials. The surface density E(r, 9,t) is Fourier 
expanded in 6: 

m — oo 

£(r,0,t) = a (r,t) + ^ a m {r,t)cos(m6 + <j) m (r,t)) (A-l) 

m— 1 

with [a m (r,t),(f>m(r,t)], the m th mode amplitude and phase 
respectively. This is practically performed with FFTs over 
a grid in 6 at various judiciously chosen radii. The relative 
power K m in a given mode m is then computed 

oo 

J a m (r y t)2irrdr 

K m {t) = , (A-2) 

J2 I a m (r,t)2nrdr 

m=0 

Cluster kinematics were determined from a snapshot of 
our benchmark 3D simulation at T=4 Myr (well into the 
saturated regime). The scale length of that simulation is on 
the order of 1 pc, whereas that of M31's nucleus is more like 3 
pc. To compare with M31 observations, we rescaled positions 
(and with them the softening length) by a factor s, and time 
by a factor si; velocities are naturally rescaled by s~2 to 
leave equations of motion in the new variables invariant. The 
adopted scaling (s ~ 2.7) was first constrained within the 
interval 2.0 < s < 3.1, by matching peaks in the rotation 
curve, then fine tuned to fit the kinematics, while adjusting 
the sky angles as well. 

While search i ng f or the optimal sky angles 
l|Peiris fc Tremainel 2003) of our rescaled snapshot, we 
took into account artifacts of the HST SITS instrument 
by convolving all the particles with a double Gaussian 
l|Bender et alj|2005h . The slit width was set to 0.1" with a 
PA of 39°. Both the fiducial, and rescaled snapshots were 
observed with the same sky angles 6 a = —59.1°, 6i = 65.3°. 
We note that larger 9i permits an almost perfect match 
with the observed LOSVD (remarkably so for 6i close to 
M31's inclination of 77°), with deteriorating rotation curve; 
while smaller 8i improved the fit to the rotation curve 
dramatically while disturbing agreement in the LOSVD. 



B THE CR INSTABILITY EXPLORED 

Here, we summarize results of a battery of experiments con- 
ducted with CR configurations that are likely to occur in 
galactic centers. 

First, we explore dynamics in highly inclined configurations, 
by considering an extreme situation: a lower resolution ( 
50xl0 3 particles ) analogue of the 3D benchmark simula- 
tion was evolved with equal mass annuli ( M r =M p =0.05 
M., M r and M p are the masses of the retrograde and pro- 
grade disks respectively) where the PP and RP disks have 
a relative inclination of 90° . This experiment was motivated 
by our concern about non-coplanar, CR rings in galactic 
centers, and whether such configuration may result in m=l 
instabilities. The configuration which was violently unstable, 
resulted in both components merging, then relaxing into a 
thick, precessing, and lopsided disk. 

Second, we explor ed the effect of a n in-spiraling cluster on 
our annular disk (|Fuiii et alj|2010h . Evolving an actual re- 
alization of a plummer sphere with Gadget-2 around the 



SMBH with enough accuracy is prohibitively expensive, thus 
the cluster was modeled as a point mass with a large soft- 
ening radius 0.1 pc. The point mass was put on a retro- 
grade Keplerian orbit around the SMBH. The massive disk 
of our 3D simulation was used as the prograde component 
(M p =0.1 M., <a>=0.92 pc), but with 25xl0 3 particles in- 
stead of 250 x 10 3 . The default values of the point mass (clus- 
ter) orbit are M c =0.01 M. (M c is the mass of the cluster), 
semi-major axis a c =0.9 pc, eccentricity e c =0, and inclina- 
tion I c =180°. The simulations for this scenario consisted of 
four suites. In each suite one parameter was varied while 
keeping the remaining parameters to their default value: 

• Suite 1. M c was varied from 0.001 to 0.1 M.: The 

prograde disk developed an m=l mode for M c > 0.01 M.. 

• Suite 2. The inclination of cluster was varied from 90° 
to 180°: In all cases the cluster and the disk ended up in 
the same plane, and with an m=l growing instability mode 
when the relative inclination was larger than 90°. 

• Suite 3. The eccentricity was varied from to 0.9: 
the low eccentricity m=l growing instability shuts off for 
e c > 0.6. 

• Suite 4. The semi-major axis was varied from 0.8 to 
10 pc: for a c > 2<a> the annulus remained axisymmetric, 
with an m=l modes growing for a c < 2<a> 
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